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I illustrate the techniques and the results of the computation of order-a s corrections to the production of 
hadrons in the collision of two virtual photons, that originate from the incoming leptons at e + e~ colliders. 



1. Introduction 

In the study of scattering phenomena in the 
limit of very large energies at fixed momentum 
transfer, one must face the fact that the two most 
promising probes of BFKL physics, namely small- 
x effects in DIS and large-rapidity-gap events in 
dijet production at hadron colliders, are inher- 
ently dependent upon long-distance physics, due 
to the presence of hadrons in the initial state; 
this poses serious difficulties on the possibility of 
cleanly extracting signals of BFKL physics from 
data. 

To overcome these difficulties a gedanken ex- 
periment has been proposed where two 
quarkonium states collide; the transverse size of 
the quarkonium is small enough to allow the 
perturbative computation of its wave function. 
Of course, since quarkonium colliders are out of 
sight, one must devise some other solutions; an 
increasing-popular one is the study of the process 



7* +7* 



hadrons, 



(1) 



at fixed photon virtualitics — qf = Q\ > 0, and 
for large center-of-mass energies squared W = 
(<7i +92) 2 ; here, we denote with the momenta of 
the photons. The virtual photons play the same 
role as the quarkonia; they are colourless, and 
their virtualities control their transverse sizes, 
which are roughly proportional to l/-y/Q|, thus 
allowing for a completely perturbative treatment. 
Notice that the virtuality of the photon is there- 
fore physically equivalent to the (squared) mass 
of the quarkonium; however, while the mass of the 
quarkonium is fixed by nature, the virtuality of 
the photon can be controlled by the experimental 
setup. 
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The easiest way to access the process of eq. ([!]) 
is through the following reaction: 



+ e~ + 7* + 7* 



hadrons: 



(2) 



namely, one considers e + e~ collisions, selecting 
those events in which the incoming leptons pro- 
duce the photons through the elementary e + e~7 
vertex, and the two photons eventually initiate 
the hard scattering that produces the hadrons. 
Data relevant to the process in eq. @ have been 
published by several experiments in the past, and 
many more are expected in the near future by 
LEP collaborations; it goes without saying that a 
linear collider would be the ideal place where to 
look at such processes. 

It is clear that the process in eq. (||) is non phys- 
ical; rather, it has to be understood as a short- 
hand notation for a subset of Feynman diagrams 
contributing to the process that is actually ob- 
served: 



e + hadrons. 



(3) 



Other contributions to the process in eq. (Q) are, 
for example, those in which the incoming e + e~ 
pair annihilates into a photon or a eventually 
producing hadrons and a lepton pair, or those in 
which one (or both) of the two photons in eq. @ 
is replaced by a Z. However, it is not difficult to 
devise a set of cuts such that the process in eq. (§) 
gives the only non-negligible contribution to the 
process in eq. (|J). One can tag both the outgo- 
ing leptons, and retain only those events in which 
the scattering angles of the leptons are small: in 
such a way, the contamination due to annihila- 
tion processes is safely negligible. Furthermore, 
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small-angle tagging also guarantees that the pho- 
ton virtualities are never too large (at LEP2, one 
typically measures Qf = O(10 GeV 2 ); therefore, 
the contribution from processes in which a pho- 
ton is replaced by a Z is also negligible. Thus, it 
is not difficult to extract the cross section of the 
process 7*7* — > hadrons from the data relevant 
to the process in eq. (||). 




while diagrams relevant to BFKL physics, such 
as d), e) and f), are expected to give 

00 

a™ ~ 1 + aj(a s Ly + 0(a s (a s L) J *), (5) 
i=i 

where L = log(W/ fi%) is a "large" logarithm, and 
all subleading logarithmic terms are indicated 
with 0(a s (a s Ly); the quantity ^ is a mass 
scale squared, of the order of the photon virtu- 
alities. By comparing eqs. and (||), it is clear 
that the latter will dominate over the former in 
the asymptotic energy region W — > 00. However, 
at current collider energies, <7 7 . 7 » is not safely 
negligible, and must be taken into proper account 
when comparing theory and data. For this rea- 
son, one usually subtracts the theoretical predic- 
tions for cr 7 « 7 * from data, and then compares the 
results obtained in this way to the predictions 
for <r^?Ji L . Unfortunately, at present only the 
leading order contribution (diagram a) in fig. |l|) 
to i7 7 * 7 * has been considered. This amounts to 
say that diagrams such as b) and c) have been 
neglected so far; these diagrams represent the 
first non-trivial QCD corrections to the process 
in eq. ([!]). We will denote these contributions 
as next-to-leading order (NLO) corrections, al- 
though effectively of leading order in a s . The 
aim of this work, which presents preliminary re- 
sults from ref. 0, is to report the computation of 
these NLO corrections. 



Figure 1. Sample of diagrams contributing to the 
production of hadrons in two-photon collisions. 



Eventually, the data are compared to the theo- 
retical predictions obtained by using BFKL equa- 
tion. Some of the diagrams contributing to the 
process in eq. ([!]) are depicted in fig. diagrams 
d), e) and f) (plus all the diagrams obtained 
by emitting more and more gluons from the t- 
channel gluon) are included in BFKL dynamics. 
In the large- W limit, diagrams with a t-channel 
quark exchange, such as a), b) and c), are ex- 
pected to give a cross section behaving as 

a rT ~ 1/W, (4) 



2. Computation of NLO corrections 

The computation of the NLO corrections to a 
hard scattering process is by now a rather stan- 
dard procedure, since algorithms exist that are 
universal (that is, process independent), and ap- 
plicable to any number of final state partons. The 
role of these algorithms is to combine in a physi- 
cally sensible way the virtual and the real contri- 
butions, that are unphysical and divergent upon 
loop and phase-space integrations. The informa- 
tion on the hard process basically enter only in 
the computation of the matrix elements. In our 
case, one needs to compute the amplitude of the 
process e + e~ — > e + e~qq at one loop, and of the 
process e + e~ — > e + e~ qqg at the tree level. Fortu- 
nately, these results are easily obtained from ex- 
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isting literature (notice that we work with mass- 
less flavours; being mainly interested in the re- 
gion of W not close to the heavy-quark thresh- 
old, our approximation is sufficiently good). As 
far as the one-loop amplitude is concerned, we 
can derive it from the one-loop amplitude rel- 
evant to the process e + e~ — > q a q a qbqb given in 
ref. Q. For the tree- level amplitude, we can 
use the results of ref. relevant to the process 
qq — > Z* Z*g —* l a l a lblb9i with a suitable crossing 
and inserting the electromagnetic coupling factors 
instead of the electroweak ones; the fact that the 
Z has an axial coupling is easily dealt with, since 
the results of ref. [Q are given in terms of helicity 
amplitudes. 





Figure 2. Decomposition of the phase space for 
the process e + e _ — » e + e~ + hadrons. 



Having the matrix elements at disposal, one 
can plug them is his preferred NLO algorithm, 
and get physical results. Our case can however 
be greatly simplified in a preliminary stage; in 
fact, the incoming and outgoing leptons actually 
do not participate in the hard scattering, that is 
initiated by the two virtual photons. Formally, 
the simplification goes through a suitable decom- 
position of the phase space, which is pictorially 
represented in fig. ||. We write the phase space of 
two leptons plus n partons (in our case, n = 2 or 
n = 3 for the one-loop or the tree-level amplitudes 
respectively) as follows 

dT{p'i,p' 2 )d^ n {qi +q 2 ;k 1 ,...,k n ), (6) 



where pi (pj) and qi = pi — p\ are the momenta 
of the incoming (outgoing) leptons and photons 
respectively, and ki are the momenta of the final 
state (coloured) partons. The first entry in d$2+n 
and e?$„ is the sum of the momenta of the incom- 
ing particles (the e + e~ pair and the two-photon 
pair respectively, as shown in fig. ||). Eq. (jfy im- 
plicitly defines dT, since we know a priori how 
to write the phase space for an arbitrary number 
of particles, given their momenta. We can thus 
immediately get 



dT 



d 3 p[ d 3 p' 2 
(2ir) 3 2p[° (2ir) 3 2p' 2 °' 



(7) 



We have therefore decomposed the original phase 
space into the product of two pieces, each of which 
is Lorentz invariant. We exploit this fact by re- 
writing eq. (0) in the center-of-mass frame of the 
incoming e + e~ pair: defining S — (pi +P2) 2 , we 
get 

dT = * dQldQ 2 2 dE 1 dE 2 dcpdip, (8) 

where ip and (p are two generic azimuthal angles, 
one of which (let's say ip to be definite) can be 
interpreted as the angle between the two outgo- 
ing leptons; are the energies of the outgoing 
leptons in the center-of-mass frame of the incom- 
ing e + e~~ pair. The strategy of the computation 
should now be clear. The hard process that we ac- 
tually deal with at NLO is that of eq. (jl]) ; thanks 
to the decomposition in eq. @, we have a 2 — > n 
phase space which is formally identical to that 
one gets as a starting point of any NLO algo- 
rithm; thus, we can safely adopt one of the ex- 
isting NLO algorithms, and study the process of 
eq. (Jl]) in the 7*7* center-of-mass frame, without 
any reference to the incoming or outgoing lep- 
tons: this amounts to a non-trivial simplification, 
since the complexity of the numerical computa- 
tions at NLO is known to grow fast with the num- 
ber of particles involved in the hard scattering. Of 
course, the information on the lepton momenta is 
entering somewhere, in particular in the matrix 
elements; to take into account this fact, we pro- 
ceed in two steps, still using fig. ^ as a guide. We 
start by generating the full kinematical configu- 
ration of the outgoing leptons, using eq. (ra). In 
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doing this, we also get the photon momenta, and 
therefore we know how to boost from the e + e~ to 
the 7*7* center-of-mass frame. Then, we boost 
the lepton momenta to the 7*7* center-of-mass 
frame, where we generate the remaining (parton) 
momenta, according to the phase space d& n ; at 
this point, we can perform all the manipulations 
required by the NLO algorithm. 

3. Results 

Following the procedure outlined in the previ- 
ous section, we constructed a user-friendly code 
capable of predicting, to NLO accuracy, any 
infrared-safe quantity constructed with up to 
three partons (plus two leptons) in the final 
state.0 The code is based upon the NLO algo- 
rithm of refs. (5||(|, and it is a suitable modifica- 
tion of one of the codes presented in ref. || . 

In order to produce phenomenological results, 
we used the following input parameters, taken 
from an analysis performed by L3 [Q; we stress 
that other physically sensible sets of parame- 
ters would lead to the same qualitative conclu- 
sions. The center-of-mass energy is fixed to \/S = 
200 GeV; the scattering angles 8i and energies 
Ei of the outgoing leptons are required to ful- 
fill 0.03 <9 t < 0.066 and E t > 40 GeV, respec- 
tively (actually, L3 use Ei > 30 GeV; we prefer 
a slightly larger value in order to have larger vir- 
tualities. The relative difference between the two 
choices is of the order of 1%). We use a two-loop 
expression for a s , with a s (M z ) — 0.1181 ||. A 
choice has also to be made for the scales entering 
the strong and electromagnetic running coupling 
constants. For the latter, we left the possibility 
open in the code to have independent scales for 
the two photon legs; it seems therefore natural 
to choose each scale equal to the virtuality of the 
corresponding leg. As far as the scale entering 
the strong coupling is concerned, we have chosen 



Ql + Ql , (HUk 



(9) 



where k T i are the transverse momenta (in 
the 7*7* center-of-mass frame) of the emitted 

2 The code can be obtained upon request by mailing to 
Stefano.Frixione@cern.ch. 



coloured partons; a discussion on this choice will 
be given in ref. ||. 

In fig. H we present our results for the total 
cross section in e + e~ collisions (that is, relevant 
to the process in eq. (g)), as a function of Y. This 
quantity is defined as follows: 



Y = log(S/S ), 5 = ^M. 



2/12/2 



where 



Vi 



1 



2Ej 



cos 2 (^/2). 



(10) 



(11) 



It is worth noticing that Y is directly related 
to the BFKL logarithm L entering eq. (||). In 
fact, for large W's, W ~ 2/i2/2<5; furthermore 
(see for example ref. ||), a sensible choice is 
Mw = c q\/ Q\Q\i with c Q a suitable constant. It 



then follows that 



L = Y - \ogc Q 



(12) 



A m-Z - 




E e+e _ = 200 GeV 
0.03<S<0.066, E„>40 GeV 



solid: NLO 



dashed: Born 



Figure 3. Cross section as a function of Y. 
LO and NLO results are both shown. 



The 



In fig. ||, the dashed histogram represents the 
LO result, while the solid histogram is the full 
NLO result; thus, the difference between the two 
histograms is the order-a s correction computed 
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in this paper. As can be seen from the figure, 
the NLO corrections are basically negligible when 
small Y's are considered. On the other hand, 
in the region of intermediate and large Y, the 
NLO corrections are definitely more important; 
at Y = 4 (Y = 6), the LO result is increased by 
a factor of about 1.2 (1.6). This is a non trivial 
piece of news; in fact, when comparing L3 data 
to LO predictions at - say - Y = 4, theory un- 
derestimates the data by a factor of about 2.5. 
This fact is interpreted as the signal of a dom- 
inant contribution coming from diagrams with 
a ^-channel gluon exchange. However, the re- 
sult given here shows that part of the discrep- 
ancy has to be associated with diagrams that do 
not have any ^-channel gluon exchange. In other 
words, the BFKL asymptotic region is probably 
still far away, and this implies that the contri- 
butions of diagrams with four quarks (that is, of 
order a§), but without i-channel gluons, is prob- 
ably not completely negligible. In any case, for 
a precise comparison between theory and data, 
NLO corrections cannot be ignored. 
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Figure 4. Dijet cross section at NLO. 



We now go back to eq. (|12[), and consider the 
role of the constant c Q . In ref. |J , this constant 
has been estimated by equating the definition of 



L with the rapidity difference between the most 
forward and backward quarks produced in the 
hard scattering. We can repeat the same exer- 
cise, using jets instead of quarks (as one always 
should, when a NLO computation is available). 
We thus computed two-jet cross sections, select- 
ing the most forward/backward jets and requir- 
ing the most (least) energetic of them to have 
transverse energy larger than 14 (10) GeV in the 
center-of-mass frame of the 7*7* pair. The differ- 
ence in their rapidities, denoted as At], is shown 
in fig. |], for various cuts on Y. It is interesting to 
notice that the contribution due to the large Y re- 
gion (dotted histogram, corresponding to Y > 6) 
basically coincides, in the large Ar/ region, with 
the cross section integrated over all Y's. We 
therefore find what we expect, namely that the 
large Y region is basically populated by events 
characterised by two hard jets well separated in 
rapidity. Thus, following ref. j9|, we can invert 
eq. (|12|) to get c Q , by substituting Y = 6 and 
identifying L with the average Arj corresponding 
to the cut Y > 6. We have 

CQ = exp(F-(Ary))~75, (13) 

which is rather close to the value of 100 estimated 
in ref. j|. The author wishes to thank Vittorio 
Del Duca for reading the manuscript. 
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